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Abstract 

We consider the electromagnetic field in a cavity with a periodically oscillating 
perfectly reflecting boundary and show that the mathematical theory of circle 
maps leads to several physical predictions. Notably, well-known results in the 
theory of circle maps (which we review briefly) imply that there are intervals 
of parameters where the waves in the cavity get concentrated in wave packets 
whose energy grows exponentially. Even if these intervals are dense for typical 
motions of the reflecting boundary, in the complement there is a positive 
measure set of parameters where the energy remains bounded. 
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I. INTRODUCTION 



In this paper, we consider the behavior of the electromagnetic field in a resonator one 
of whose walls is at rest and the other moving periodically. The main point we want to 
make is that several results in the mathematical literature of circle maps immediately yield 
physically important conclusions. 

The problem at hand is mathematically equivalent to the study of the motion of vibrating 
strings with a periodically moving boundary [|T],0, or the classical electromagnetic field in a 
periodically pulsating cavity fHH]. It is connected with the vacuum quantum effects in such 
region |5|||. The problem is also of practical importance, e.g., for the formation of short 
laser pulses @. 

The goal of this paper is to show that the problem of a classical wave with a periodically 
moving boundary can be easily reformulated in terms of the study of long term behavior 
of circle maps and, therefore, that many well known results in this theory lead to physical 
predictions. In particular, we give proofs of several results obtained numerically by Cole 
and Schieve |3J] and others. Extensions of this approach, that will be discussed elsewhere, 
allow to reach conclusions for some quasi-periodic motions of small amplitude or possibly 
for non-homogeneous media. 

In the case of more than one spatial dimension, the analogous problem |J is much more 
complicated, so the predictions are not as clear as in the one- dimensional case and we will 
not discuss them further. 

We emphasize that the mathematical theory presented is completely rigorous and, hence, 
the physical predictions made are general for the assumptions stated. 

There are other intriguing relations for which we have no conceptual explanation. We 
observe that a calculation of Fulling and Davies || leads to the conclusion that the energy 
density radiated by a moving mirror is equal to the Schwarzian derivative of the motion of 
the mirror (for details see Sec. [IV E|) . This Schwarzian derivative is a differential operator 



frequently used in the theory of one dimensional dynamical systems and particularly in the 
theory of circle maps. 

The plan of the exposition is the following. In Sec. |I| we show how the physical problem 
can be formulated in terms of circle maps. Sec. |T| contains a brief exposition of the necessary 
facts from the theory of circle maps, in Sec. [IV] these facts are applied to the problem at 
hand and illustrated numerically, and in the conclusion we discuss the advantages of our 
approach. 

II. PHYSICAL SETTING 

A. Description of the system 

We consider a one-dimensional optical resonator consisting of two parallel perfectly re- 
flecting mirrors. For simplicity of notation, we will consider only the situation in which one 
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of them is at rest at the origin of the x axis while the other one is moving periodically with 
period T. The case where the two mirrors are moving periodically with a common period 
can be treated in a similar manner. We assume that the resonator is empty, so that the 
speed of the electromagnetic waves in it is equal to the speed of light, c. The speed of the 
moving mirror cannot exceed c. 

We note that the experimental situation does not necessarily require that there is a 
physically moving mirror. One experimental possibility - among others - would be to 
have a material that is a good conductor or not depending on whether a magnetic field 
of sufficient intensity is applied to it, and then have a magnetic field applied to it in a 
changing region. This induces reflecting boundaries that are moving with time. Note that 
the boundaries of this region could move even faster than c, hence the study of mirrors 
moving at a speed comparable to c is not unphysical (even if, presumably, one would also 
have to discuss corrections to the boundary conditions depending on the details of the 
experimental realizations). 

We shall use dimensionless time t and length i connected with the physical (i.e., dimen- 
sional) time tph y s and length £ phys by t := t phys /T, t := £ phys /(cT). 

Let the coordinate of the moving mirror be x = a(t), where a is a C k function (k = 
1, . . . , oo, uS) satisfying the conditions 

a{t)>0, \a'{t)\<l, a(t + l) = a(t). (2.1) 

The meaning of the first condition is that the cavity does not collapse, the second one means 
that the speed of the moving mirror cannot exceed the speed of light, and the third one is 
that the mirror's motion is periodic of period 1. An example which we will use for numerical 
illustrations is 

a(f) = | + P sin27rt (l/3|<^> 0<|/3|<|). (2.2) 

Since there are no charges and no currents, we impose the gauge conditions Aq = 0, 
V • A = on the 4-potential = (A , A) and obtain that A satisfies the homogeneous 
wave equation. We consider plane waves traveling in x-direction, so that without loss of 
generality, we assume that A(t,x) = A(t,x)e y , and obtain that A(t,x) must satisfy the 
homogeneous (1 + l)-dimensional wave equation, 

A u (t,x)-A xx (t,x) = 0, (2.3) 

in the domain £ := {(t, x) G R 2 \ t < t, < x < a(t)}. It will also need to satisfy some 
boundary conditions that will be specified in Sec. [II B| , and appropriate initial conditions, 



A(t ,x) =ip 1 (x), A t (t ,x) =ip 2 {x) . (2.4) 

Before discussing the boundary conditions and the method of solving the boundary- value 
problem in the domain E, let us discuss the way of solving ( |2.3| ) in the absence of spatial 
boundaries, i.e., in the domain {t Q < t, x G R}. It is well-known that in this case, the 
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solution of the problem (|2.3|) , (|2.4|) at some particular space-time point (t, x) can be written 

as 



A(t,x) = ^ 



Xn 



I Xp 



(2.5) 



where ± (t — t ), and ^ and are functions of one variable that are se- 

lected to match the initial conditions (|2.4j ). The explicit expressions for ^ ± follow from 
the d'Alembert's formula (see, e.g. 



0) 

Ms) ± 



Ms') da' 



(2.6) 



where £ is an arbitrary constant (the same for and $ ). 

The representation Q2.5Q has a simple geometrical meaning: the value of A(t, x) is a 



superposition of two functions, ^~(x$ ) and \l/ + (x ( j"), the former being constant along the 
lines {x — t = const}, and the latter being constant along {x + t = const}. The disturbances 
at a space-time point (T, X) propagate in the space-time diagram along the lines {x — t = 
X — T} and {x + t = X + T} emanating from this point (in more physical terms, this 
corresponds to two rays moving to the right and to the left at unit speed); these lines are 
called characteristics, and the method of solving ( |2.3| ), (|2.4j ) by using the representation 
(2.5) is called the method of characteristics (see, e.g., |10|,|TTf ). 



B. Method of characteristics, boundary condition at the moving mirror, and Doppler 

shift at reflection 

To obtain the boundary conditions at the stationary mirror, we note that the electric 
field, i.e., the temporal derivative of the vector potential, must vanish at this mirror, which 
yields the following "perfect reflection" boundary condition: 

AM) = 0. (2.7) 

The boundary condition at the moving mirror can be easily obtained by performing a 
Lorenz transformation from the laboratory frame to the inertial frame comoving with the 
moving mirror at some particular moment t. The temporal and spatial coordinates in the 
laboratory frame, t and x, are related to the ones in the instantaneously comoving inertial 
frame, t' and x', by 

t — to = t! cosh ( + x' sinh ( 

(2.8) 

x — a(to) = t! sinh ( + t' cosh ( , 

where tanh£ = a'{t). In the comoving frame, the boundary condition is A t >(0, 0) = 0, which, 
together with (|2.8j), yields 

A t (t, a{t)) + a' (t) A x {t, a(t)) = . (2.9) 
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The method of characteristics developed in ( |2.5| ) and ( |2.6| ) for situations with no bound- 
aries can be adapted to provide rather explicit solutions for systems in spatially bounded 
space-time domains satisfying (^3J) at the boundaries (see, e.g., (T^J). 

The prescription is the following. The solution of the boundary value problem ( |2.3| ), 
(|2.4D, (|2.7|), ( ^.9[ ) in the domain £ is a superposition of two functions that are constant on 
the straight pieces of the characteristics and change their sign at each reflection. To find 
A(t,x), one has to consider the two characteristics, 7" and 7" 1 ", passing through (t, x), and 
propagate them backwards in time (according to the rule that, upon reaching a mirror, they 
change direction of propagation) until they reach the line {time = to} a t the points (to, Xq ) 
and (to, a^o ), resp. - see Fig. [|. Then A(t,x) is given by 

A(t,x) = (-l^-tf-fo) + (-l)"+¥+0n+) , (2.10) 

where are the number of reflections of 7 T on the way back from (t,x) to (to,^). In 
Sec. [11 Q we will give explicit formulae for x^ and A(t, x) in terms in circle maps. 



Indeed, because the solution (|2.10|) is the sum of two functions constant along the straight 
pieces of the characteristics, the wave equation is satisfied in the interior. Also, the initial 
conditions are easily verified because for t — to small, Xq and Xq are close to x [see ( p,17|) 1. 



To check that this prescription also satisfies the boundary conditions, we need another 
argument. Consider the space-time diagram of the reflection of the field between two in- 
finitesimally close characteristics reflected by the moving mirror at time 9, shown in Fig. |2|. 
The world line of the mirror is denoted by m, the angle 5 between it and the time direction 
is connected with the mirror's velocity at reflection by tan 6 = a' (9). The Doppler factor 
at reflection, D(9), is defined as the ratio of the spatial distances A and A' between the 
characteristics before and after reflection: 

n/M A /7T A l-tan5 1 - a'(9) 

D(9 := — = tan 5 ) = = H" • (2.11) 

A' V4 J l + tan5 1 + a'{9) V ; 

Thus, the absolute values of the temporal and spatial derivatives of the field increase by a 
factor of D(9) after reflection. This implies that if in the space-time domain between the 
two characteristics, the values of the corresponding derivatives of the field before reflection 
are denoted by A t and A x , then after reflection they will become —D(9)A t and D(9)A X , resp. 
Hence, in the space-time domain of the overlap the derivatives of the field will be 

A t (9,a(9)) — A t - D(9)A t 

(2.12) 

A x (9,a(9)) = A x + D(9)A x . 
Now, we will show that the modified method of characteristics is consistent with the bound- 



ary condition (|2.9| ). We note that A t = —A x , which simply means that before reflection the 
rays are moving to the right at unit speed. If we multiply the second equation of (|2.12|) 
by a!(9) = [which follows from ([2.1 1|) ] and add it to the first, we obtain exactly the 



boundary condition (fO) 
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The same prescription gives a solution of the Dirichlet problem A(t, 0) = A(t, a{t)) = 0. 
Similar methods can be developed for other boundary conditions. 

We note that the method of characteristics also yields information in the important 
case when the medium is inhomogeneous and perhaps time dependent. This is a physically 
natural problem since in many applications we have cavities filled with optically active media 
whose characteristics are changed by external perturbations. In this case, the method of 
characteristics does not yield an exact solution as above but rather, it is the main ingredient 
of an iterative procedure JTT[. Physically, what happens is that in inhomogeneous media, 



the waves change shape while propagating in contrast with the propagation without change 
in shape in homogeneous media (|2.5|). We plan to come to this problem in a near future. 



C. Using circle maps to solve the boundary- value problem fl2.3|) , (|2.4j ), ( |2.7[) , (2.9) in 

the domain £ 

We now reformulate the method of characteristics into a problem of circle maps. 

We consider a particular characteristic and denote by {r n } the times at which it reaches 
the stationary mirror and {0 n } the times at which it reaches the oscillating one; let ... < 
T n < d n < r n+ i < 6 n+ i < . . .. Note that, with this notation, 

r n = 0„-a(0 n ) = (Id-a)(0 n ) 

(2.13) 

r n +i = n + a(0„) = (ld+a)(0 n ) . 

Therefore 

r n+1 = (Id+a) o (Id-a) _1 (r n ) =: F(r n ) 

(2.14) 

6 n+1 = (Id -a)" 1 o (ld+a)(0„) =: G(6 n ) . 

We refer to F and G as the time advance maps. They allow to compute the time of reflection 
on one side in terms of the time of the previous reflection on the same side. The conditions 
( [2.1D on the range of a and a 1 guarantee that (Id —a) is invertible and that F and G are C k 
(by the implicit function theorem). 

When the function a is 1-periodic, F and G satisfy 

F(t + l) = F(t) + l, G{t + l) = G{t) + l. (2.15) 

These relations mean that F(t) and G(t) depend only on the fractional part of t. In physical 
terms, we characterize a reflection of a ray by the phase of the oscillating mirror when the 
impact takes place, i.e., by the time of reflection modulo 1; if we know the phase at one 
reflection, we can compute the phase at the next impact. Mathematically, this means that 
F and G can be regarded as lifts of maps from S 1 = R/Z to S* 1 (see Sec. P~TT|) . 

We want to argue that the study of the dynamics of the circle maps ( |2.14j ) leads to 
important conclusions for the physical problem, which we will take up after we collect some 
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information about the mathematical theory of circle maps. In particular, many results in 
the mathematical literature are directly relevant for physical applications. This is natural 
because the long term behavior of the solution can be obtained by repeated application of 
the time advance maps [see Q2.18| )]. 



We call attention to the fact that 

G = (Id +a) -1 o F o (Id +a) = (Id -a)" 1 o F o (Id -a) , (2.16) 

so that 

G n = (Id +a)' 1 oF n o (Id +o) = (Id -ay 1 oF n o (Id -a) . 

In dynamical systems theory this is usually described as saying that the maps F and G 
are "conjugate" (see Sec. |III C|) . In our situation, this comes from the fact that F and G 
are physically equivalent descriptions of the relative phase of different successive reflections: 
F advances the r variables while G advances the #'s, and the #'s are related to the r's by 



Now, we use circle maps to derive an explicit formula for the solution of the boundary- 
value problem (2.3), (2.4), (|2.7f) , (|2.9| ) in the domain E. Let us trace back in time the 



characteristics 7" and 7 + coming "from the past" to the space-time point (t, x) - see Fig. [l|. 
Let 6^ := (a =F Id) _1 (t — x) be the last moments the characteristics 7 ± are reflected by 
the moving mirror, and let 9^ k := G~ k (9o). After iV + , resp. AL, reflections on the way 
backwards in time (out of which n+, resp. n_, are from the moving mirror), the characteristic 
7 + , resp. 7", crosses the line {time = to}- The spatial coordinate of the intersection of 7 ± 
and {time = t } can be easily seen to be 

4 = h(9± n± , t ) := \(ld-a)(6%J - to\ • (2.17) 

Thus, the formula for the vector potential is 

A(t, x) = (-1)*-*- o h (G~ n - o (a + ld)-\t - x), t Q ) 

+ (-!)*+*+ o h (G- n+ o (a - Id)- 1 ^ - x), t ) . (2.18) 



D. Energy of the field 



The method of characteristics gives a very illuminating picture of the mechanism of the 
change of the field energy, 



E(t) 



a(t) 



T 00 (t,x) 



■Sir Jo 



a(t) 



A t (t,x) 2 + A x (t,x) 2 dx, 



(2.19) 



due to the distortion of the wave at reflection from the moving mirror. Indeed, consider the 
change of the energy of a very narrow wave packet at reflection from the moving mirror at 
time 6. Since at reflection the temporal and the spatial distances decrease by a factor of 
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D(8), \A t \ and \A X \ will increase by a factor of D(9). Therefore, the integrand of the energy 
integral will increase D{9) 2 times, while the support of the integrand (i.e., the spatial width 
of the wave packet at time t) will shrink by a factor of D(8). Hence, the energy of the wave 
packet after reflection will be D(9) times greater than its energy before reflection. 

In the general case, one can use (J2.18]) and obtain the energy of the system at time t. 



For the sake of simplicity, we will give the formula only under the assumption that at time t 
all the rays are going to the right, i.e., assuming that the vector potential is of the form 
A(t,x) = (—1) n -^~(xq). Let us introduce the "local Doppler factor" 



D(t , x ; t) :-- 



d 

—h(9 n _, t ) 



l + a'(9 ) 1 



(2.20) 



It has the physical meaning of the ratio of the frequencies of the incident wave and the wave 
at time t [cf. (|2~Tl| )1. Note that D(t Q , Xq ; t) is equal to the derivative of G n ~ multiplied by 
a factor which is bounded and bounded away from independently of [due to the fact 
that \a'(t)\ < 1]. From ( |2.18|) and (|2.17|) we obtain that the square of D(t , Xq ; t) is the 
ratio of the energy density T 00 (t,x) and the initial energy density, T 00 (t , Xq ): 

T 00 (t,x) = 2 |(^-)'(x -)| 2 D(t , x ; tf = T 00 (t , ar") D(t Q , x "; t) 2 . 
On the other hand, D(to, x^; t) is connected with the Jacobian of the change of coordinates 



Xn 



x by 



dx 




dx 


dxo 




dx 



D(t , x ; t) 



-i 



Hence, the energy of the system at time t is 

r a(t) 



E(t) = / r 00 (t , xq) D{t , Xq- t) dxo 
Jo 



(2.21) 



Note that since the local Doppler factor squared is the ratio of the energy densities at 
two consecutive reflection points, then it satisfies the following multiplicative property. Let 
(tx,Xi), (#2, ^2 )> • • • > be space-time points on the characteristic connecting (t ,XQ ) 

and (t, x), such that at all of them the rays are going to the right, and let to < t\ < . . . < 
tk <t. Then 

D(t , Xq \ t) = D(t , x ; ti) D(ti, x\\ t 2 ) ■■■ D(t k -i, x k _ x \ t k ) D(t k , x k ; t) . 

As can be seen from (|2.20|) , these multiplicative properties are closely related to the chain 
rule for diffeomorphisms, 



(G n )'{9) = G'(G n -\9)) G'{G n ~ 2 {9)) ■ ■ ■ G'(9) 



(2.22) 



The mathematical theory of dynamical systems contains many results about derivatives 
of highly iterated maps as above ( |2.22| ). In Sec. [IV C| we will be able to translate some of 
them into asymptotic properties of the field energy. 
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A simple and intuitively clear formula for the rate of change of the field energy can be 
obtained by using ( [2.19] ), Q2.3|), ( |2.7| ), (|2.9|), and integrating by parts: 

1 Mt) 1 
= T~ / + A - A ^ dx + 7rCL'{t)A x {t, a{t)) 2 

47T JO 07T 

= —A t (t,a(t))A x (t,a(t)) + ^-a'{t)A x {t,a{t)) 2 

1 * „ 



57T 



: a'(t)A x (i,a(t)) 2 . 



From this, we realize that the force experienced by the wall is ^A x (t, a{t)) 2 . 

E. The inverse problem: determining the mirror's motion given the circle map 

It is important to know whether the notion of a "typical" G is the same as the notion of 
a "typical" a or a "typical" F (in the mathematical literature people speak about "generic" 
maps, and in physical literature about "universal" maps). We do not know the answer to 
this question, and here we will give some arguments showing that the answer is not obvious. 
In this paper we will not use "generic" or "universal" . Rather we will make explicit the non- 
degeneracy assumptions so that they can be checked in the concrete examples. In Sec. [IV D 



we will show that some universal properties for families of circle maps do not apply to G 
constructed according to fl2.14| ) with a(t) = a + eb(t). 

While the function a can be expressed in terms of F as a = F ~ ld o ( F ~^ j , the relation 
between G and a is much harder to invert. We should have 

a{9) + a{G{8)) = G{9) , (2.23) 

where G{6) := G{6) — 6, so for any n, 

a(9) = G{6) - G{G{6)) + ■■■ + {-l) n G{G n {6)) + (-l) n+1 a(G n+1 (9)) . 

Hence, if G 2k (6o) = 6q (modi), a necessary condition for the existence of a is that 

2fe-l 

y £(-l) i G(G i (0o)) = 0. (2.24) 



i=0 



An example of a G where the above condition is not satisfied can be readily constructed. We 
note that if a map fails to satisfy ( |2.24 ) and if (G 2k )'(6 ) ^ 1, then all small perturbations 



will also fail to satisfy (|2.24j) . Thus, there are whole neighborhoods of maps that cannot be 
realized as G for a moving mirror. 

On the other hand, given very simple G's, it is easy to construct infinitely many a's 
that satisfy ( |2.23|) and that therefore lead to the same G. For example, for G(8) = 9 + |, 



( p.23j ) amounts to a{6 + |) + a(9) = |. If we prescribe a for 9 in [0, |], then this equation 
determines a on [|, 1] (the only care needs to be exercised so that the two determinations of 
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a match at 9 = 5). A similar construction works when G permutes several intervals - if we 
prescribe a(9) in an interval J, ( |2.23| ) determines a{9) in G(I). 



In the case when G is conjugate to an irrational rotation, G = h~ l o R a o h, then ( |2.23| ) 
is equivalent to 

a o h^ 1 o R + a o /tT 1 = h~ l o i? a — /tT 1 . 



Then ao/i 1 can be determined using Fourier analysis, setting h 1 (9) = 9 + J2kL-oo %e 2mke , 
a o br x (d) = 9 + ET=-oo ^ke 2mke , which leads to 

^brita + 1 ^ k= ^Mka _ ~ (2 25) 

Let us assume that \ka — n — || > const l&l -1 ' for some v > 1 (a condition of this type is 
called a Diophantine condition - see definition |III.4|) , and that h~ x has r derivatives (which 



implies that its Fourier coefficients % satisfy \r k \ < const |fc|~ r ). Then if r > v + 2, then 
the coefficients ipk define a smooth function (for more details see, e.g., Sec. XIII.4]). Of 
course, once we know aoh" 1 , then, since h~ x depends only on G and is therefore determined, 
we can obtain a. 

In summary, there are maps G that do not come from any a at all, come from infinitely 
many a's, or come from one and only one a. The maps F can always be obtained from one 
and only one a. 



III. MAPS OF THE CIRCLE 



In this section we recall briefly some facts from the theory of the dynamics of the orien- 
tation preserving homeomorphisms (OPHs) and diffeomorphisms (OPDs) of the circle S 1 , 
following closely the book of Katok and Hasselblatt [14], Ch. 11, 12]; see also and |fl3| . 



This is a very rich theory and we will only recall the facts that we will need in the physical 
application. 

We shall identify S* 1 with the quotient R/Z and use the universal covering projection 



7T 



: R -> S 1 = R/Z : x h+ n(x) := x (mod 1) . 



Another way of thinking about S* 1 is identifying it with the unit circle in C, and using the 
universal covering projection x 1— > e 2mx . 

Let / : S 1 —>■ S 1 be an OPH and F : R — > R be its lift to R, i.e., a map satisfying 
/ o 7r = 71 o F. The fact that / is an OPH implies that F(x + 1) = F(x) + 1 for each x G R, 
which is equivalent to saying that F — Id is 1-periodic. The lift F of / is unique up to an 
additive integer constant. If a point x G S 1 is g-periodic, i.e., f q (x) = x, then F q (x) = x + p 
for some p G N. 
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A. Rotation number 



A very important number to associate to a map of the circle is its rotation number, 
introduced by Poincare. It is a measure of the average amount of rotation of a point along 
an orbit. 

Definition III.l Let f : S 1 — > S 1 be an orientation preserving homeomorphism and F : 
R — > R a lift of f . Define 

T (F) := Km ~ x , r(/ ) := To (F) (mod 1) (3.1) 

and call r(f) a rotation number of f. 

It was proven by Poincare that the limit in ( |3.1| ) exists and is independent of x. Hence, 
r(/) is well defined. 

The rotation number is a very important tool in classifying the possible types of behavior 
of the iterates of the OPHs of S 1 . The simplest example of an OPH of S 1 is the rotation 
by a on S 1 = R/Z, r a : x i— > x + a (modi) (corresponding to a rotation by 2na radians 
on S 1 thought of as the unit circle in C). The map R a : x i— > x + a is a lift of r a , and 
T ( r a) = « (modi). In the case of r a there are two possibilities: 

(a) If r(r a ) = p/g G Q, then R^ q (x) = x + p for each x e R, so every point in S 1 is 
g-periodic for r p / q . If p and g are relatively prime, q is the minimal period. 

(b) If r{r a ) ^ Q, then r a has no periodic points; every point in S 1 has a dense orbit. 
Thus, the a- and u;-limit sets of any point x G 5* 1 are the whole S 1 , which is usually 
described as saying that S 1 is a minimal set for r a . [Recall that a(x) is the set of the 
points at which the orbit of x accumulates in the past, and u(x) those points where it 
accumulates in the future.] 



B. Types of orbits of OPHs of the circle 



To classify the possible orbits of OPHs of the circle, we need the following definition (for 
the particular case / : S 1 — ► S 1 ). 

Definition III. 2 (a) On orbit O of f is called homoclinic to an invariant set T G S l \0 
if a(x) = u(x) = T for any x G O. 

(b) An orbit O of f is said to be heteroclinic to two disjoint invariant sets T\ and T 2 if O 
is disjoint from each of them and a(x) = T\, u(x) = T 2 for any x G O. 

With this definition, the possible types of orbits of circle OPHs were classified by Poincare 



T6| as follows (for a modern pedagogical treatment see, e.g., [0, Sec. 11.2]) 
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(1) For r(f) = p/q G Q, all orbits of / are of the following types: 

(a) a periodic orbit with the same period as the rotation r p / q and ordered in the same 
way as an orbit of r p / q ; 

(b) an orbit homoclinic to the periodic orbit if there is only one periodic orbit; 

(c) an orbit heteroclinic to two different periodic orbits if there are two or more 
periodic orbits. 

(2) When r(f) ^ Q, the possible types of orbits are: 

(a) an orbit dense in S 1 that is ordered in the same way as an orbit of r r (/) (as are 
the two following cases); 

(b) an orbit dense in a Cantor set; 

(c) an orbit homoclinic to a Cantor set. 

We also note that in cases 2(b) and 2(c), the Cantor set that has a dense orbit is unique 
and can be obtained as the set of accumulation points of any orbit. 



C. Poincare and Denjoy theorems 



Because of the simplicity of the rotations it is natural to ask whether a particular OPH of 
S l is equivalent in some sense to a rotation. To state the results, we give a precise definition 
of "equivalence" and the important concept of topological transitivity. 

Definition III.3 (a) Let f : M —> M and g : N N be C m maps, m > 0. The maps f 
and g are topologically conjugate if there exists a homeomorphism h : M —> N such 
that f = h^ 1 o g oh. 

(b) The map g is a topological factor of f (or f is semiconjugate to g) if there exists a 
surjective continuous map h : M — > iV such that h o / = g o h; the map h is called a 
semiconjugacy. 

(c) A map f : M — *> M is topologically transitive provided the orbit, {f k {x)}k&z, of some 
point x is dense in M . 

The meaning of the conjugacy is that g becomes / under a change of variables, so that 
from the point of coordinate independent physical quantities, / and g are equivalent. The 
meaning of the semiconjugacy is that, embedded in the dynamics of /, we can find the 
dynamics of g. 



The following theorem of Poincare |T6| was chronologically the first theorem classifying 
circle maps. 

Theorem III.l (Poincare Classification Theorem) Let f : S 1 — > S 1 be an OPH with 
irrational rotation number. Then: 
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(a) if f is topologically transitive, then f is topologically conjugate to the rotation ?>(/)/ 

(b) if f is not topologically transitive, then there exists a non-invertible continuous mono- 
tone map h : S 1 — > S 1 such that ho f = r T (/) o h; in other words, f is semiconjugate 
to the rotation r T /f) ■ 

If we restrict ourselves to considering not OPHs, but OPDs of the circle, we can say 
more about the conjugacy problem. An important result in this direction was the theorem 
of D enjoy Jlj 



Theorem III. 2 (Denjoy Theorem) AC 1 OPD of S 1 with irrational rotation number and 
derivative of bounded variation is topologically transitive and hence (according to Poincare 
theorem) topologically conjugate to a rotation. In particular, every C 2 OPD f : S 1 — > S 1 is 
topologically conjugate to r T m . 



We note that this condition is extremely sharp. For every e > there are C 2 e maps 



with irrational rotation number and semiconjugate but not conjugate to a rotation (see |13|, 
Sec. X.3.19]). 

D. Smoothness of the conjugacy 

So far we have discussed only the conditions for existence of a conjugacy h to a rotation, 
requiring h to be only a homeomorphism. Can anything more be said about the differentia- 
bility properties of h in the case of smooth or analytic maps of the circle? As we will see 
later, this is a physically important question since physical quantities such as energy density 
depend on the smoothness of the conjugacy. To answer this question precisely, we need two 
definitions. 

Definition III. 4 A number p is called Diophantine of type (K, v) (or simply of type v) for 
K > and v > 1, if p — 2 > K \q\^ 1 ^ u for all - G Q. The number p is called Diophantine 
if it is Diophantine for some K > and v > 1. A number which is not Diophantine is called 
a Liouville number. 

It can be proved that for K — > 0, the set of all Diophantine numbers of type (K, v) has 
Lebesgue measure as close to full as desired. 

Definition III. 5 A function f is said to be C m ~ s where m > 1 is an integer and 5 & (0, 1), 

if it is C m ~ x and its (m — l)st derivative is (1 — 5)-Holder continuous, i.e., 



D m - l f{x) - D m - l f{y) 



< const \x — y] 1 5 
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The first theorem answering the question about the smoothness of the conjugacy was the 
theorem of Arnold ||18|| . He proved that if the analytic map / : S 1 — > S 1 is sufficiently close 
(in the sup-norm) to a rotation and r(f) is Diophantine of type v > 1, then / is analytically 
conjugate to the rotation r r (/), i.e., there exists an analytic function h : S 1 — > S 1 such that 
h o f — r T (j) o /i. The iterative technique applied by Arnold was fruitfully used later in the 
proof of the celebrated Kolmogorov- Arnold- Moser (KAM) theorem - see, e.g., jn|. Arnold's 



result was extended to the case of finite differentiability by Moser ||20|| . In such a case, the 
Diophantine exponent v has to be related to the number of derivatives one assumes for the 
map. 

Arnold's theorem is local, i.e., it is important that / is close to a rotation. Arnold conjec- 
tured that any analytic map with a rotation number in a set of full measure is analytically 
conjugate to a rotation. Herman |TJ| proved that there exists a set A C [0, 1] of full Lebesgue 
measure such that if / G C k for 3 < k < u> and r(f) G A, then the conjugacy is C k ~ 2 ~ £ 
for any e > 0. The set A is characterized in terms of the growth of the partial quotients 
of the continued fraction expansions of its members; all numbers in A are Diophantine of 
order v for any v > 1. His result was improved by Yoccoz who showed that if / G C k , 



3 < k < lj, r(f) is a Diophantine number of type v > 1, and k > 2v — 1, then there exists 
a C k ~ u ~ £ conjugacy h between / and r T m for any e > 0, and by several others. The best 
result on smooth conjugacy we know of, is the following version of Herman's theorem as 
extended by Katznelson and Ornstein p2[. 



Theorem III. 3 (Herman, Katznelson and Ornstein) Assume that f is a C k circle 
OPD whose rotation number is Diophantine of order v, and k > v + 1 . Then the homeo- 
morphism h which conjugates f with the rotation r T (n is of class C k ~ v ~ e for any e > 0. 

There are examples of C 2 ~ e maps with a Diophantine rotation number arbitrarily close 
to a rotation and not conjugated by an absolutely continuous function to a rotation - see, 

e.g., pa. 



E. Devil's staircase, frequency locking, Arnold's tongues 

Let {f a }aeA be a one-parameter family of circle OPHs such that f a (x) is increasing in 
a for every x. Then the function a i— > r(f a ) is non-decreasing. (Since the maps are only 
defined modulo an integer and so is the rotation number, some care needs to be taken to 
define increasing and non-decreasing when some of the objects we are considering change 
integer parts. What is meant precisely is that if one takes the numbers with their integer 



parts, they can be made increasing or non- decreasing. This is done in detail in [13, Sec. 11.1], 
and we will dispense with making it explicit since it does not lead to confusion.) 

For such a family the following fact holds: if r(f a ) ^ Q, then a t— > r(f a ) is strictly 
increasing locally at a; on the other hand, if f a has rational rotation number and the 
periodic point is attracting or repelling (i.e., there is a neighborhood of the point that gets 
mapped into itself by forwards or backwards iteration), then a >— > r(f a ) is locally constant 
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at this particular value of a, i.e., for all a' sufficiently close to a, r(/ a /) = r(/ a ). The local 
constancy of the function a i— > r(/ a ) is known as frequency (phase, mode) locking. Note 
that, since the rotation number is continuous, when it indeed changes, it has to go through 
rational numbers. The described phenomenon suggests the following definition. 

Definition III. 6 A monotone continuous function ip : [0, 1] — > R is called a devil's staircase 
if there exists a family {/g}^ G = of disjoint open subintervals of [0, 1] with dense union such 
that ip takes constant values on these subintervals. (We call attention to the fact that the 
complement of the intervals in which the function is constant can be of positive measure.) 

The devil's staircase is said to be complete if the union of all intervals 1^ has a full 
Lebesgue measure. 

A very common way of phase locking for differentiable mappings arises when the map 
we consider has a periodic point and that the derivative of the return map at the periodic 
point is not equal to 1. By the implicit function theorem, such a periodic orbit persists, and 
the existence of a periodic orbit implies that the rotation number is locally constant. At the 
end of the phase locking interval the map has derivative one and experiences a saddle-node 
(tangent) bifurcation. 

We note that, unless certain combinations of derivatives vanish (see, e.g., ||24||), the 
saddle-node bifurcation happens in a universal way. That is, there are analytic changes of 
variables sending one into another. This leads to quantitative predictions. For example, the 
Lyapunov exponents of a periodic orbit should behave as a square root of the distance of 
the parameter to the edge of the phase locking interval. 

Of course, other things can happen in special cases: the fixed point may be attractive but 
only neutrally so, there may be an interval of fixed points, the family may be such that there 
are no frequency locking intervals (e.g., the rotation). Nevertheless, all these conditions are 
exceptional and can be excluded in concrete examples by explicit calculations. (For example, 
if the family of maps is analytic but not a root of the identity, it is impossible to have an 
interval of fixed points.) 

In the example we will consider, we will not perform a complete proof that a devil's 
staircase occurs, but rather we will present numerical evidence. In particular, the square 
root behavior of the Lyapunov exponent with the distance to the edge of the phase locking 
interval seems to be verified. 

Let us now consider two-parameter families of OPDs of the circle, {4> a ^}, depending 
smoothly on a and (3. Assume that when (3 = 0, the maps of the family are rotations by a, 
i.e., Qj o = T a . We will call (3 the nonlinearity parameter. Assume also that dtpa^/da > 0. 
An example of this type is the family studied by Arnold 



Va,/3 '■ S 1 — > S 1 : x i— > 7] a ,p{x) : = x + a + (3 sin 2-kx (mod 1) , (3.2) 

where a e [0,1), (3 e (0,1/2tt). 

The rotation number r is a continuous map in the uniform topology, and <p a ^ is a con- 
tinuous function of a and (3, so the function (a, (3) i— > r(0 Q , /3 ) =: Tp(a) depends continuously 
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on a and (3. The map Tp is non-decreasing; for (3 > 0, Tp is locally constant at each a for 
which Tp(a) is rational and strictly increasing if Tp(a) is irrational. Thus, Tp is a devil's 
staircase. 

Since Tp is strictly increasing for irrational values of Tp(a), the set I u := {(a, (3) \ Tp(a) = 
u} for an irrational v G [0, 1] is a graph of a continuous function. For a rational u, I u has 
a non-empty interior and is bounded by two continuous curves. The wedges between these 
two curves are often referred to as Arnold's tongues. 

The fact that r(0 Q , j o) = T~(r a ) = a implies that for (3 = 0, the set of a's for which there 
is frequency locking coincides with the rational numbers between and 1, so its Lebesgue 
measure is zero. When (3 > 0, its Lebesgue measure is positive. The width of the Arnold's 



tongues for small (3 for the Arnold's map ( |3.2| ) is investigated, e.g., in p3 |. Much of this 
analysis carries out for more general functions such as the ones we encounter in the problem 
of the periodically pulsating resonator. 

The total Lebesgue measure of the frequency locking intervals, m({Tp l (v) \ v G Q fl 
[0, 1]}), becomes equal to 1 when the family of circle maps consists of maps with a horizontal 
point (so that the map, even if having a continuous inverse, fails to have a differentiable 



one) - see |£6[27| f° r numerical results and [28] for analytical proof. With the Arnold's map 



rj a ^p this happens when (3 = l/2n. In our case this happens when the mirror goes at one 
instant at the speed of light. 

We note also that the numerical papers P^j26|j2^j30|] contain not only conjectures about 
the measure of the phase locking intervals but, perhaps more importantly, conjectures about 
scaling relations that hold "universally" . In particular, the dimension of the set of parameters 
not covered by the phase locking intervals should be the same for all non-degenerate families. 
These universality conjectures are supported not only by numerical evidence but also by a 
renormalization group picture - see, e.g., |^TJ and the references therein. These universality 
predictions have been verified in several physical contexts. Notably in turbulence by Glazier 



and Libchaber 32 



As we will see in Sec. [IV D| , we do not expect that the families obtained in ( |2.14| ) for 
mirrors oscillating with different amplitudes belong to the same universality class as typical 
mappings, but they should have universality properties that are easy to figure out from those 
of the above references. 



F. Distribution of orbits 

For the physical problem at hand it is also important to know how the iterates of the 



circle map x i— > g{x) := G(x)(modl) are distributed. As we shall see in lemma [IV. 1| , if 
the iterates of g are well distributed (in an appropriate sense), the energy of the field in the 
resonator does not build up. The distribution of an orbit is conveniently formalized by using 
the concept of invariant measures. We recall that a measure fi on X is invariant under the 
measurable map / : X — > X if fi(f~ l (A)) = fi(A) for each measurable set A. 

Given a point ifS 1 , the frequency of visit of the orbit of x to / C S 1 can be defined by 
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m r #{E|0<i<nand f (x) e 1} , . 

/,,(/) := hm . (3.3) 



It is easy to check that if for every interval /, the limit ( |3.3|) exists, it defines an invariant 
measure describing the frequency of visit of the orbit of x. Therefore, if there are orbits 
which have asymptotic frequencies of visit, we can find invariant measures. 

A trivial example of the existence of such measures is when x is periodic. In such 
a case, the measure \i x is a sum of Dirac delta functions concentrated on the periodic 
orbit. The measure of an interval is proportional to the number of points in the orbit it 
contains. We also note that it is easy to construct two-dimensional systems [e.g., (r, 6>) i— > 
(1 + 0.1(r — 1), 9 + (r — l) 2 sin(#) 2 )] for which the limits like the one in (|3.3| ) do not exist 
except for measures concentrated on the fixed points, so that even the existence of such 
equidistributed orbits is not obvious. 

There are also relations going in the opposite direction - if invariant measures exist, 
they imply the existence of well distributed orbits. We recall that the Krylov-Bogolyubov 
theorem Jl4|, Thm. 4.1.1] asserts that any continuous map on a compact metrizable space has 



an invariant probability measure. Moreover, the Birkhoff ergodic theorem [jy|, Thm. 4.1.2] 
implies that given any invariant measure /x, the set of points for which \i x as in (|3.3|) does 
not exist has measure zero. 

Certain measures have the property that \i x — fi for //-almost all points. These measures 
are called ergodic. (There are several equivalent definitions of ergodicity and this is one of 
them.) From the physical point of view, we note that a measure is ergodic if all the points in 
the measure are distributed according to it. For maps of the circle, there are several criteria 
that allow to conclude that a map is ergodic. 

For rotations of the circle with an irrational rotation number we recall the classical 
Kronecker-Weyl equidistribution theorem [|T4| , Thm. 4.2.1] which shows that any irrational 



rotation is uniquely ergodic, i.e., has only one invariant measure - the Lebesgue measure m. 
(Such uniquely ergodic maps are, obviously, ergodic because, by Birkhoff ergodic theorem, 
the limiting distribution has to exist almost everywhere, but, since there is only one invariant 
measure, all these invariant distributions have to agree with the original measure.) Thus, 
the iterates of any x G S 1 under an irrational rotation are uniformly distributed on the 
circle. 

For general non-linear circle OPDs the situation may be quite different. As an example, 
consider Arnold's map r] aj p ( |3.2| ). If it is conjugate to an irrational rotation by h, i.e., 
Va,p = ° r T{r) afj ) ° h, then there is a unique invariant probability measure /i defined for 
each measurable set A by fJ>(A) := m(h(A)). This implies that if I is an interval in S 1 , then 
the frequency with which a point x visits / is equal to //(/). 

On the other hand, if r(r] a ^) =p/g £ Q, then all orbits are periodic or asymptotic to 
periodic. Thus, the only possible invariant measure is concentrated at the periodic points 
and therefore singular, if the periodic points are isolated. Let us now assume that a is very 
close to Tg (p/q), but does not belong to it. Then i] a ^ has no periodic orbits, but still there 
exists a point x which is "almost periodic", i.e., the orbits linger for an extremely long time 
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near the points x, r] a)/3 (x), ■ ■ • , tj^~}(x). So that, even if the invariant measure is absolutely 
continuous, one expects that it is nevertheless quite peaked around the periodic orbit - see 
Fig. The behavior of such maps is described quantitatively by the "intermittency theory" 



The continuity properties of the measures of the circle are not so easy to ascertain. 
Nevertheless, there are certain results that are easy to establish. 

Of course, in the case that we have a rational rotation number and isolated periodic 
orbits, some of them attracting and some of them repelling, the only possible invariant 
measures are measures concentrated in the periodic orbits. 

For the irrational rotation number case, the Kronecker-Weyl theorem implies that all the 
maps with an irrational rotation number - since they are semi-conjugate to a rotation by 
Poincare theorem - are uniquely ergodic. In the situations where Herman's theorem applies, 
this measure will have a smooth density since it is the push-forward of Lebesgue measure 
by a smooth diffeomorphism. 

We also recall that by Banach-Alaoglu theorem and the Riesz representation theorem, 
the set of Borel probability measures is compact when we give it the topology of //„ — > 
fi /i n (^4) - ► M^) f° r a ^ Borel measurable sets A. (This convergence is called weak-* 

convergence by functional analysts and convergence in probability by probabilists.) 

Lemma III.l If X* is a parameter value for which f\* admits only one invariant measure 
fix*, given invariant measures for f\ i , with \ — > X* , then fii converges in the weak-* 
sense to 

Note that we are not assuming that f\ t are uniquely ergodic. In particular, the lemma 
says that in the set of uniquely ergodic maps, the map that a parameter associates the 
invariant measure is continuous if we give the measures the topology of weak-* convergence. 

Proof. Let /iA ifc be a convergent subsequence. The limit should be an invariant measure for 
fx* . Hence, it should be . It is an easy point set topology lemma that for functions taking 
values in a compact metrizable space, if all subsequences converge to the same point, then 
this point is a limit. The space of measures with weak-* topology is metrizable because by 
Riesz representation theorem is the dual of the space of continuous functions with sup-norm, 
which is metrizable. 



We also point out that as a corollary of KAM theory |18| we can obtain that for non 



degenerate families, if we consider the parameter values for which the rotation number 
is Diophantine with uniform constants, the measures are different iable jointly on x and 
in the parameter. (For the differentiability in the parameter, we need to use Whitney 
differentiability or, equivalently, declare that there is a family of densities differentiable 
both in x and in A that agrees with the densities for these values of A.) 

On the other hand, we point out that there are situations where the invariant measure 
is not unique (e.g., a rational rotation or a map with more than one periodic orbit). In 
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such cases, it is not difficult to approximate them by maps in such a way that the invari- 
ant measure is discontinuous in the weak-* topology as a function of the parameter. The 
discontinuity of the measures with respect to parameters, as we shall see, has the physical 
interpretation that, by changing the oscillation parameters by arbitrarily small amounts, we 
can go from unbounded growth in the energy to the energy remaining bounded. 

IV. APPLICATION TO THE RESONATOR PROBLEM 

Now we return to the problem of a one-dimensional optical resonator with a periodically 
moving wall to discuss the physical implications of circle maps theory, and illustrate with 
numerical results in an example. 



A. Circle maps in the resonator problem 

If we take a(t) to depend on two parameters, a and (3, as in (|2.2j), then, as we saw in 



Sec. |II C| , the time between the consecutive reflections at the mirrors can be described in 
terms of the functions F Q /3 and G a ^ defined by (|2.14j) . These maps are lifts of circle maps 
that we will denote by and g a ^. The restriction on the range of (3 in ( |2.2| ) implies that 
f a) p and Qol,(3 are analytic circle OPDs. Therefore, we can apply the results about the types 
of orbits of OPHs of S 1 , Poincare and Denjoy theorems, as well as the smooth conjugacy 
results and the facts about the distribution of orbits. 

In an application where the motion of the mirror [i.e., a(t)] is given, one needs to compute 
-F Qj/3 and G a ^ Q2.14D , which cannot be expressed explicitly from a(t) but they require only to 
solve one variable implicit equation. In the numerical computations we used the subroutine 
zeroin |34| to solve implicit equations. If y = F a> p(t) and z = G a ,p(t), then for a(t) given 



by (p.2|), y and z are given implicitly by 



- y + t + a + 2(3 sin[vr(?/ + t)] = 
- z + t + a + (3 [sin(27rt) + sin(27r^)] = 



Given t, we can find y, z applying zeroin. 



B. Rotation number, phase locking 

In this section, our goal is to translate the mathematical predictions from the theory of 
circle maps into physical predictions for the resonator problem. 

The theory of circle maps guarantees that the measure of the frequency locking intervals 
for g a> p is small when (3 is small and becomes 1 when (3 = l/2n. The theory also guarantees 
for analytic maps that, unless a power of the map is the identity, the frequency locking 
intervals are non-trivial. For the example that we have at hand, it is very easy to verify 
that this does not happen and, therefore, we can predict that there will be frequency locking 
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intervals and that as the amplitude of the oscillations of the moving mirror increases so that 
the maximum speed of the moving mirror reaches the speed of light, the devil's staircase 
becomes complete. Fig. ^ shows a part of the complete devil's staircase - the situation which 
happens when the maps g a ^ and f a ^ lose their invertibility, i.e., for (3 = l/2n. 

We also recall that the theory of circle maps makes predictions about what happens for 
non-degenerate phase locking intervals. Namely, for parameters inside the phase locking 
interval the map has a periodic fixed point and the Lyapunov exponent is smaller than 0, 
while at the edges of the phase locking interval the map experiences a non-degenerate saddle- 



node bifurcation - provided that certain combinations of the derivatives do not vanish [24 



We note that for parameters for which the map is in non-degenerate frequency locking, 
i.e., r{g a ^) = p/q and the attractive periodic point of period q has a negative Lyapunov 
exponent, {G™^ /3 (x)}^L will converge exponentially to the fixed point for all x in a certain 
interval, according to the results about the types of orbits of circle maps (Sec. |111B| ). The 
whole circle can be divided into such intervals and a finite number of periodic points. There- 
fore, the graph of G^ q g, and hence of g™^, will look - up to errors exponentially small in n - 
like a piecewise-constant function with values (up to integers) in the fixed points of g q a @ - see 
Fig. [|. The fact that certain functions tend to piecewise-constant functions for large values 
of the argument (which follows from what we found about G n £a for large n) was observed 
numerically for particular motions of the mirror in H|J. In physical terms, this means that 
the rays will be getting closer and closer together, so with the time the wave packets will 
become narrower and narrower and more and more sharply peaked. The number of wave 
packets is equal to q. The number of reflections from the moving mirror per unit time will 
tend to the inverse of the rotation number. In the next section we discuss how this yields 
an increase of the field energy which happens exponentially fast on time. 

The fact that for r(g a ^) G Q the rays approach periodic orbits, is also interesting from 
a quantum mechanical point of view due to the relation between the periodic orbits in a 
classical system and the energy levels of the corresponding quantum system, given by the 
Gutzwiller's trace formula (see, e.g., p5| ). 

We also note that we expect that slightly away from the edges of a phase locking interval, 
the invariant density will be sharply peaked around the points in which it was concentrated 
in the phase locking intervals. This is described by the "intermittency theory" . 



To observe numerically in our example what happens when a enters or leaves a fre- 
quency locking interval, we set Np{y) := {a G [0, 1) | r(g al3 ) = u}. Fig. || represents the 
Radon-Nikodym derivative d/z/dm of the invariant probability measure \i with respect to the 
Lebesgue measure m [i.e., of the density of the invariant measures, which, as we saw in ( |3.3|) , 
is the frequency of visit of the iterates]. The figure shows d/x/dm for a close to the left end 
of iV ai (i/6). When a approaches (from the left) the left end of iV 01 (l/6), dfi/dm becomes 
sharply peaked at some points, and when a enters the frequency locking interval, the invari- 
ant measure becomes singular (g Q ,o.i undergoes tangent bifurcation at a = 0.253977...). 
All seems to be consistent with the conjecture that all the frequency-locking intervals in the 
family (away of (3 = 0) are non-degenerate, i.e., that at the boundaries of the phase locking 
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intervals the map satisfies the hypothesis of the saddle-node bifurcation theorem. 



C. Doppler shift 

One of the most interesting parts of the applications of circle map theory is the ease with 
which we can describe the effect on the energy after repeated reflections. 



Recall that in Sec. [II D| , we found the time dependence of the field energy under the 
assumption that at time t all rays are going to the right. This assumption is not very 
restrictive in the case of a rational rotation number since, as we found in Sec. [LVB , the field 



develops wave packets that become narrower with time, so (|2.20|) and ( |2.21| ) hold for the 



asymptotic behavior of the energy. Note that ( 2.20 ) expresses the Doppler shift factor in 
terms of the derivatives of the map G. This gives a very close relation between the dynamics 
and the behavior of the wave packets. 

Proposition IV. 1 Let a and j3 be such that r(g a ^) = p/q, and that the map G := G a> p has 
a stable periodic orbit <d q = {9\, . . . , 9 q } such that (G q )'(9\) < 1. Assume that the initial 
electromagnetic field in the cavity is not zero at some space-time point for which the phase 
of the first reflection from the moving mirror is in the basin of attraction of Q q . 

Then the energy of the field in the resonator will be asymptotically increasing at an 
exponential rate: 



Proof. First notice that the number of reflections from the moving mirror per unit time 
reaches a well defined limit (one and the same for all rays) - the inverse of the rotation 
number. Secondly, as was discussed in Sec. ||, at reflection from the moving mirror at 
phase 9, a wave packet becomes narrower by a factor of D{9) (|2.11| ), which leads to a 
D{9) times increase in its energy. Asymptotically, the phases at reflection will approach the 
stable periodic orbit Q q = {6\, . . . , 9 q } of g a ^. The Doppler factors at reflection will tend 
correspondingly to {D(9i), . . . , D{9 q )} (|2.11| ). Hence, in time p each ray will undergo q 



reflections from the moving mirror, the total Doppler shift factor along the periodic orbit 
Q q being 

On the other hand, the definition of the map G as the advance in the time between 
successive reflections from the moving mirror yields 9i = G l ~ 1 {9{). The chain rule applied 
to the explicit expression fl2.14|) for G yields 



9-1 9-1 l + a'(9- 

{G q ~ i y{9i) = n g\9 3 ) = n 

, = 1 7=1 L ~ a W-+ 
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which gives the following expression for D(Q g ) [cf. ( |2.20|) 1: 

D(B q ) := [(G^M)] ~ l = \ Z ^M^)\6 q ) . (4.2) 

1 + a'{9 q ) L J 1 + a'(9 q ) 

Hence, the energy density grows by a factor of D(Q q ) 2 . Since after q reflections the wave 
packet is concentrated in a length D(O q ) times smaller, the total energy grows by a factor 
of D(Q q ) in p units of time, which implies fl4.1|) . 



The quantities (G n )'(9) that appear in (|4.2| ) have been studied intensively in dynamical 
systems theory since they control the growth of infinitesimal perturbations of trajectories. 
Similarly, they are factors that multiply the invariant densities when they get transported, 
as we will see in ( j4.3|) . 

We found numerically the total Doppler factors D(Q q ) for some particular choices of 
the parameters. In Fig. ||, log 10 -D(O 6 ) is shown for different values of f3 and for a G 
Np(l/6). Obviously, the maximum value of D(Qq) depends strongly on f3, becoming infinite 
for (3 = l/2n and some a G -/Vi/ 27r (l/6). For smaller values of (3, the Doppler factor is 
much smaller. Moreover, the width of the frequency locking intervals for small (3 is small, 
so the probability of hitting a frequency locking interval with arbitrarily chosen a and (3 is 
small. [The likelihood of frequency locking for the Arnold's map ( p.2| ) is studied numerically 
in P.] 

In the case when Herman's theorem apply, the derivatives of G n are bounded indepen- 
dently of n, which causes the energy of the system to be bounded for all times, which is 
proved in the following proposition. 

Proposition IV. 2 If G a ^ is such that it satisfies the hypothesis of Herman's theorem, then 
the energy density remains bounded for all times. 

Proof. In such a case G a ^ = hT 1 o R o h with h differentiable and R a rotation by r(g a ^). 
Therefore G™ ^ = h^ 1 o R n o h and 

(G n a J(e) = (h-')'(R" o h{9)) (R n )'(h(9)) h'(9) = (h- 1 )'^ o h{9)) h'(9) 

because (R n )' = 1. The two factors in the right-hand side of the above equation are bounded 
uniformly in 9 and n. Thus, the "local Doppler factors" ( |2.20| ) will be bounded, which implies 
the boundedness of the energy (|2.21|). 



There is an interesting connection between the invariant densities of the system and the 
growth of the electromagnetic energy density. 

Recall that if a density \i is invariant, fi(G(9)) = ji{9)/G'{9). Hence, if the density 
H never vanishes, G'{9) = ^i{9)/^{G{9)) and, therefore, (G i )'(9) = fi(9) / /i(G< (9)) . Let us 
assume that there is only one characteristic passing through the space-time point (t, x), and 
this characteristic is going to the right. Then, using the notations of Sec. [11 Q we can write 
the energy density at (t, x) as [cf. Q2.20Q 1 
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r 00 (t,x) 



l-q'(0Z n _) fi(G n -(ez n _)) 
1 + a'(9o) ^(9Z n _) 



T 00 (t ,x ). 



(4.3) 



In the general case [with two characteristics through (x, t)], one can use fl2.18|) and ([2.1 9|) 
to prove the following result: 



Lemma IV. 1 If a system has an invariant density /i which is bounded away from zero, then 
the electromagnetic energy density of a C 1 initial data is smaller than Cfi 2 for all times. 

In the cases that Herman's theorem applies, there is an invariant density bounded away 
from zero (and also bounded). Hence, we conclude that there are values of the amplitude of 
mirror's oscillations for which the energy density of the field remains bounded. This set is 
typically a Cantor set interspersed with values for which the energy increases exponentially. 

Some other results about the behavior of the energy with respect to time and parameters 
are obtained in 0. 

We call attention to the fact that [18| contains examples of analytic maps whose rotation 
numbers are very closely approximated by rationals and that are arbitrarily close to a rota- 
tion such that they preserve no invariant density and, therefore, are not smoothly conjugate 
to a rotation. 

It is also known that for all rotation numbers one can construct C 2 ~ e maps arbitrarily 
close to rotations with this rotation number and such that they do not preserve any invariant 
It is a testament to the ubiquity of these maps that these questions were 
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motivated and found applications in the theory of classification of C* algebras. 



D. The behavior for small amplitude and universality 



We note that, even if all the motions of the mirror lead to a circle map as in ( |2.15| ), it 
does not seem clear to us that all the maps of the circle can appear as F, G for a certain a. 
This makes it impossible to conclude that the theory of generic circle maps applies directly 
to obtain conclusions for a generic motion of the mirror. Of course, all the conclusions of the 
general theory that apply to all maps of the circle apply to our case. Those conclusions that 
require non-degeneracy assumptions will need that we verify the assumptions. Nevertheless, 
the very developed mathematical theory of generic or universal circle maps cannot be applied 
without caution to maps that appear as the result of generic or universal oscillations of the 
mirror. 

One aspect that we have found makes a big difference with the generic theory is the 
situation where the mirror oscillates with small amplitude, i.e., a £ (t) = a + eb(t) with b 
a periodic function of zero average and period 1, and e < 1. The first parameter, a, is 
the average length of the resonator, while e = is called the "nonlinearity parameter" for 
obvious reasons. If we denote by F^ £ and Ga, £ the corresponding 2-parametric families of 
maps of the circle constructed according to ( 2.14| ), then we have, for three times differentiable 
families, 
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F Se (t) =t + 2a + 2eb(t + a) + 2e 2 b'(t + a)b(t + a) + 0(e 3 ) , 

(4.4) 

G a ,M =t + 2a + e[b{t) + b{t + 2a)} + e 2 b'{t + 2a)[b{t) + 6(t + 2a)] + 0{e 3 ) . 

Note that the term of order e has vanishing average. As we will immediately show, this 
property causes that some well known generic properties of families of circle mappings do 



not hold for families of maps constructed as in Q2.14D . 

Indeed, if we consider the expressions for small amplitude developed in (4.4), we can 
write the maps as 

H e (t) =t + 2a + eHi{t) + e 2 H 2 (t) + 0(e 3 ) . 

Since the conclusions of the theory of circle maps are independent of the coordinate system 
chosen, it is natural to try to choose a coordinate system where these expressions are as 
simple as possible. Hence, we choose h £ (t) := t + er](t), a perturbation of the identity, and 
consider h~ x o H £ o h £ , which is just H £ in another system of coordinates, related to the 
original one by h £ . Then, up to terms of order e 3 , we have 

h; 1 oH £ o h £ {t) =t + 2a + e [rj(t) - rj(t + 2a) + Hi{t)\ 

(4.5) 

+e 2 {q'it + 2a)r)(t + 2a) - jf{t + 2a) [r)(t) + H^t)] + H[(t)r)(t) + H 2 (t)} . 

We would like to choose i] in such a way that the e term is not present. Note that since 
J rj(t + 2a) dt = J r)(t) dt, this is impossible unless / H^t) dt = 0. When / H x (t) dt = 0, R x 



is smooth and 2a is Diophantine, a well-known result (see, e.g., [|T3j, Sec. XIII. 4]) shows that 
in such a case we can obtain one r\ satisfying 

r){t) - r](t + 2a) + H x (t) = (4.6) 

and Tj = 0. [Such rj is conventionally obtained by using Fourier coefficients. Note that in 
Fourier coefficients, ( |4.6| ) amounts to rjk(e 2mk2a — 1) = {Hi) k . If Hi is smooth, the Fourier 
coefficients decrease fast and if 2a is Diophantine, then ( e 27r * fc2a _ i)^ 1 does not grow too 
fast. For more details we refer to the reference above.] 

Since for the functions Fa t£ and Ga, e the term of order e has a zero average, we can 
transform these functions into lifts of rotations plus 0(e 2 ). This implies, in particular, that 
their rotation number is t(F Sj£ ) = r(G S|e ) = 2a + 0(e 2 ). One could wonder if it would be 
possible to continue the process and eliminate also to order e 2 . 

If we look at the e 2 terms in ( |4.5|) , we see that r)'(i)r)(t) = 0, and, when r] is chosen as in 



(EH) 



rf(t + 2a)[r](t) + H x (t)\ = r}'(t + 2a)r](t + 2a) , 

which also has average zero. Therefore, a necessary condition for the e 2 term in h~ l oH £ oh £ {t) 
to be zero is H[(t)rj(t) + H 2 (t) = 0. 

For the F Si£ in ( |4.4|) we see that F 2 has zero average. Nevertheless, the term F[{t)r]{t) 
does not in general have average zero as can be seen in examples. Hence, we see that the 
rotation number indeed changes by an order which is 0(e 2 ) and not higher in general. This 
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property is not generic for families of circle maps starting with a rotation 2a and it puts 
them outside of the universality classes considered in etc., since the correspondence 

between rotation numbers and parameters is not the same. 

the space of 



According to the geometric picture of renormalization developed in BT 



circle maps is divided into slices of rational rotation numbers, which are - in appropriate 
sense - parallel. In that language - in which we think of families of circle maps as curves in 
the space of mappings - the families of advance maps Fa, e and G a , e (for fixed a) have second 
order tangency to the foliation of rational rotation numbers rather than being transversal. 
Hence, the scaling predicted by universality theory should be true for e 2 in place of e. We 
have not verified this prediction, but we expect to come back to it soon. 



E. Schwarzian derivative in the problem of moving mirrors 



Fulling and Davies §| calculated the energy-momentum tensor in the two-dimensional 
quantum field theory of a massless scalar field influenced by the motion of a perfectly re- 
flecting mirror (see also [0). They obtained that the "renormalized" vacuum expectation 
value of the energy density radiated by a moving mirror into initially empty space is 



T 



00, 



u 
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u 



'F"(u) 



where u = t — x, and F is related to the law of the motion of the mirror, x = a(t), by ( p,14[ ). 
The right-hand side of this equation is nothing but (up to a constant factor) the Schwarzian 
derivative of F - a differential operator which naturally appears in complex analysis, e.g., 
it is invariant under a fractional linear transformation; vanishing Schwarzian derivative of 
a function is the necessary and sufficient condition that the function is fractional linear 
transformation, etc. More interestingly, the Schwarzian derivative has been used as an 
important tool in the proof of several important theorems in the theory of circle maps - see, 
e.g., pT| , |37| . In the light of the connection between the solutions of the wave equation in a 
periodically pulsating domain and the theory of circle maps it is not impossible that this is 
not just a coincidence. 



V. CONCLUSION 

Using the method of characteristics for solving the wave equation, we reformulated the 
problem of studying the electromagnetic field in a resonator with a periodically oscillating 
wall into the language of circle maps. Then we used some results of the theory of circle maps 
in order to make predictions about the long time behavior of the field. We found that many 
results in the theory of circle maps have a directly observable physical meaning. Notably, 
for a typical family of mirror motions we expect that the electromagnetic energy grows 
exponentially fast in a dense set of intervals in the parameters. Nevertheless, it remains 
bounded for all times for a Cantor set of parameters that has positive measure. 
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There are several advantages of the approach presented here. First, it allows us to 
understand better the time evolution of the electromagnetic field in the resonator and the 
mechanism of the change in the field energy. Second, the predictions are based on the 
general theory of circle maps so they are valid for any periodic motion of the mirror; let 
us also emphasize that our method is non-perturbative. Last, but not least, for a given 
motion of the mirror, one can easily make certain predictions about the behavior of the field 
by simply calculating the rotation number of the corresponding circle map, and without 
solving any partial differential equations. 
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FIG. 1. Finding A(t, x) by the method of characteristics. 
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FIG. 2. Reflection by the moving mirror. 



30 



0.33 





1 1 1 1 


1 

1/3 

' 4/13 

— 3/10 

« 2/7 

— 3/11 




1/4 

• 3/13 

- ' 2/9 

1/5 

2/11 




1/6 






1/7 







0.34 0.36 0.38 0.4 0.42 

a 

FIG. 3. A part of the graph of T(g a l / 27r ) vs. a. 
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FIG. 4. Development of the piecewise-constant structure of go?2545,o.i (the rotation number of 
50.2545,0.1 is 1/6)- Graphs of 5^2545,0.1 are plotted for n = 1 (dotted line), n = 5 (dashed line), 
n = 10 (long dashed line), n = 100 (solid line). 
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FIG. 5. Density of the invariant measures for [3 = 0.1 and a = 0.253 (dashed line), a = 0.2539 
(solid line), and a = 0.253975 (dotted line). 
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FIG. 6. A log-linear graph of the total Doppler factor for in the phase locking interval 
of rotation number 1/6 for different (3. The insert [linear-linear graph of D(Qq) vs. a — a c ] calls 
attention to the square-root behavior at edges; a c is the value of a at the left end of iVo.i4(l/6). 
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